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Interference effects in tunnelling of ’’cat” wave packet states 
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We analyse tunnelling of a single particle, whose initial state is given by a superposition of spatially 
separated wave packet modes. It is shown that ’’pile up” of different components the scatterer may 
change the tunnelling probabilities, making such states a convenient tool for probing the barrier’s 
scattering times. Interference effects arising in resonance tunnelling in are studied in detail. The 
analysis allows us to gain further insight into the origin of interference effects in scattering of several 
identical particles. 
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I. INTRODUCTION 


Recently, there has been considerable interest in the in¬ 
terference effects accompanying scattering of several non¬ 
interacting identical particles m-m- In the celebrated 
Hong-Ou-Mandel setup (see, for example, 0 : 01 ): such 
effects are observed if the particles, incident from oppo¬ 
site sides, coincide in the scatterer. Recently it was shown 
that a different kind of interference effects may arise in 
the case, where identical particles, incident from the same 
side, and detained in the scatterer, coincide there due to 
a kind of ’’pile up” effect [1], [7]. As a result, in reso¬ 
nance tunnelling, the transmission probability was shown 
to oscillate as a function of the temporal delay between 
the arriving particles, provided two or more metastable 
states in the barrier can be accessed [6]. 

In this paper, we discuss a closely related case, where the 
state of a single particle, incident on a barrier, consists 
of several spatially separated wave packet modes. Such 
exotic ’cat’ state can be created, for example, by splitting 
the original wave packet into parts, which experience dif¬ 
ferent time delays before being recombined [S] or, in the 
case of cold atoms, by using techniques similar to those 
described in [9]. We will show that the ’’pile up” of the 
modes, caused by a delay in the barrier region, can cause 
observable changes the tunnelling probability. One pur¬ 
pose of this paper is to analyse the use of such systems 
as an alternative tool for probing the barrier’s scattering 
times. Its other purpose is to use this analysis in order to 
gain further insight into the nature of interference effects 
in scattering of several identical particles. 

The rest of the paper is organised as follows. In Sect. I we 
consider transmission of a multi-component initial state. 
In Sect. II we analyse its transmission across a rectan¬ 
gular barrier. In Sect. Ill we study the case of resonance 
tunnelling, and the interference patterns occurring in the 
transmission probability. In Sect. IV we discuss the in¬ 
terference mechanism, and its similarity with the case of 
several identical particles. Sect. V consider transmis¬ 
sion of a mixed ’cat’ state, and Sect. VI contains our 
conclusions. 
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FIG. 1. (Color online) Schematic diagram showing a ’cat’ 
state, consisting of N non-overlapping components, incident 
on a resonance barrier supporting two metastable states 


II. A MULTI-COMPONENT INITIAL STATE 

Consider, in one dimension, a particle whose wave 
function is given by a superposition of N wave packets, 
delayed in time relative to each other {h = I), ipn{x), 

N 

^'o(x,f) = ^■!/;„(a:,f), (1) 

n—1 

ipn{x,t) = j Ari(p)eyiY>[ipx-iE{p){t + tn)]dp, 

where 0 = ti < t 2 < ... <%. For a particle of mass 
p, e.g., for cold atoms [5] or photons in a waveguide m, 
the energy is quadratic in the momentum, E(p) = p^/2p. 
For massless particles, e.g., free photons, or electrons in 
graphene m, this relation is linear, E{p) = cp. The 
constituent wave packets may, or may not overlap, and 
for the normalisation constant A in Q we have 

K = = ( 2 ) 

mn 

dpA*^{p)An{p) excp[iE{p)Tmn] = E d-mri') 

m,n 

where Tmn = tm — tn- The particle is incident on a finite 
width potential barrier with a transmission amplitude 
T{p) (see Fig.I) and, as t > oo, its transmitted part 
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takes the form 

<) = ^ V'J(a;, t) = x (3) 

n 

/ T{p)An{p) exp[ipx - iE{p){t + tn)]dp, 

n ^ 

From Eq.([^ the transmission probability we have 

(4) 

m,n m,n 

where 

Tmn = (ClV-n) = (5) 

j dp\T{p)f A*^{p)Anip) exp[iE{p)Tmn], 

is the matrix of the overlaps between the transmitted 
components 

If the delays are large, \tm — tn\ —>- 00 , rapid oscillations 
of the exponentials in Eqs.([^ and ([^, make off-diagonal 
elements of the overlap matrix vanish, Imn = o-ndmn and 
d^mn ~ J \d^{p) l^n(p)| dp X If SO, all 

components of 'kg are transmitted independently, and 
we have 

P^ = Y,Wn = Pl^. (6) 

n 

As m |t)j . w^e are interested m the case w^here lynn — 
cindmn, and Tmn ^ Wn^mn- This would indicate that 
initially non-overlapping components of 410 ’’pile up” in 
the barrier region, and the interference between them 
affects the outcome of the tunnelling process. A deviation 
of from P^ii may, therefore, serve as a crude indicator 
that a scattered particle spends in the barrier a duration 
comparable to at least some of \tm —tn\- Next we apply 
this test to the case of a rectangular barrier. 
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FIG. 2. (Color online) Tunnelling probabilities for the rect¬ 
angular barrier Q vs. v/to, to = l/2^(6 — a)^, with 
2p,V{b — a)^ = 4. Shown are the cases of A = 2 (solid) 
and N = 5 (dashed) Gaussian components ( refyla) with 
Po{b—a) = 1.41, and a/{b—a) = 4.47. To the right the vertical 
dashed line the overlap between the components, \Ej\, is 
less than 0.005. Also shown in the inset are |A(p)p (solid) and 
|r(p)|^|A(p)p (dashed), both renormalised to unit heights, as 
well as |T(p)| (thick solid). 

■0^ is shifted towards higher p’s, but not modified suf¬ 
ficiently to prevent integrals in Eq. © from being de¬ 
stroyed by the oscillations (see inset in Fig. 2). There is, 
therefore, no evidence that different components of the 
wave function, well separated initially, may be delayed, 
and eventually ’’meet” in the barrier region. The trans¬ 
mission probability for an 2- and 5-component states is 
shown in Fig.2. As expected, as soon as the overlap be¬ 
tween different ipn vanishes, different modes in Eq.([^ 
tunnel independently, and we have P^ = Pf^d- 
sence of the said pile up effect is consistent with the origi¬ 
nal McColl’s suggestion that there is no appreciable delay 
in tunnelling across a rectangular barrier [12]. It is also 
consistent with the finding of Ref. [Bj , where tunnelling of 
two identical particles was studied in a similar context. 



III. A RECTANGULAR BARRIER 


IV. RESONANCE TUNNELLING 


For a rectangular barrier, V{x) = V for a < x << 
b, and 0 otherwise, the transmission coefficient in the 
tunnelling regime E{p) < U is given by the well known 
expression 

|r(p)P = I/{1 + sinh2[g(6 - a)]/4U(U - E{p)}, (7) 

where q{p) = [2p,{V — E)]^/^ for a massive particle. 
Henceforth, we will consider Gaussian wave packets with 
identical momentum distributions, separated by equal 
time delays, 

AniP) = Amip) = A{p), tn+l-tn=T, (8) 

\Mp)\‘^ = (27r)"^/Vexp[-(p-po)^crV2]. 

In the deep tunnelling regime, |T(p)p ^ exp[—2^(6— 
a)] rapidly grows as E increases, but contains no sharp 
features. As a result, the momentum distribution of each 


The situation is different in resonance tunnelling across 
a symmetrical barrier. The transmission coefficient of 
such barrier typically exhibits well separated sharp nar¬ 
row peaks which, in the Breit-Wigner approximation, 
have a Lorentzian shapes. 




r2 
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(pV2/r - E-y + r2 


(9) 


Thus, Ej = E{pj) gives the position of the j-th reso¬ 
nance peak, and Fj is its width, and the transmitted mo¬ 
mentum distribution, |r(p)p|A(p)p|, can be made much 
narrower than the incident one, |A(p)p|. Approximat¬ 
ing both |A(p)p and dpE constant for E{p) « E”, and 
evaluating the remaining integrals in (§, then yields 

Tmn = E exp(—F jIto — n|r) exp[iE'J(m — n)T] (10) 

3 
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FIG. 3. (Color online) Probability to tunnel with one reso¬ 
nance level accessible to a particle in 2-component initial state 
Q vs. the time lag r = t 2 — ti, for poa = 6, EI/Eq = 1, 
Fi/ifo = 0.014 and Eo = E{po). A horizontal dashed line 
marks in Eq.(|^; ipi and '02 can be considered non¬ 

overlapping to the right of the vertical dashed line. The inset 
shows |A(p)p and \T{p)\‘^ 


where Cj = 27rTj\A{pj)\‘^/dpE{pj). Eq. (10) shows that 
Tmn(T) oscillates with the internal frequencies of the res¬ 
onance barrier, luj = Ej, j = 0,1,2,.... In particular, 
for N = 2, and just one resonance state at we have 
the interference correction = P^{t) — 

given by (t EE T12) 


5P^{t) = 


2pLT:Ti 

dpE{pl 


■|^(pDPexp(-rir)cos(E[r). (11) 


For a narrow resonance, SP'^{t) in Eq.( 11) oscillates with 


a frequency E^, persists for the delays at which 0i and 
02 no longer overlap, and finally vanishes for r exceeding 
l/Ei, as illustrated in Fig. 3a. This is an agreement with 
the broadly accepted view that, in resonance tunnelling, 
a particle spends approximately a duration of order of 
the life time of the metastable state supported by the 
barrier(see also [H]). 

With only two resonances accessible to the incident par¬ 
ticle, for Cl ^ C 2 = C, and |'r(ri — r 2 )| << 1, we find 


5P"^[t) « 2C'exp(—E it) cos(dwr) cos(ajr). 


( 12 ) 


where uJ = {E{ + ^0/2, and dw = {E{ + E^)/2. If 
two resonance levels are close to each other, uj » Soj, 
damped rapid oscillations of SP{t) are modulated with 
a much lower frequency 6uj (see Fig. 3b). We recall that 
for two identical particles, quantum statistical correction 
to transmission probability was found to oscillate with 
the frequency 26u! [ 3 , which suggests certain similarity 
between two effects. We will return to discuss it in Sec¬ 
tion VI. 

Finally, for an initial state ([^, containing N identical 
modes, we have 


FIG. 4. (Color online) Same as Fig.3, but for two resonance 
levels, poa = 9, E^/Eq = 0.9, Fi/Fio = 0.032, E[/Eo = 1,1, 
and Vi/Eo = 0.038. Also shown by a thick solid line is the 
envelope 2C'exp(—Fir) cos(5ajr) in Eq.(12l. 



FIG. 5. (Color online) Same as Fig. 4, but for five identical 
Gaussian components, N = 5, and tn — C-i = r. The inset 
shows the sequences of peaks [cf. ( |14[ )], contributed by each 
of the two resonances. 


where the function F/v(F’",r) is given by 

= (14) 

exp[f£J {N - l)r] - exp(-z£jT) 
[exp(-z£jT) - I ]2 

and we have introduced complex energies = FJ — iVj 
to shorten the notations. For N » 1, it is sufficient 
to retain only the first term in the curly brackets, which 
yields 


Fn{t,EAT) = -CN 


cos{E''t) — exp(— Ft) 
cos{E^t) — cosh(rr) 


Thus, for Tr << 1, Fjv(r, F’’, T)) has sharp peaks at 
r = 27rfc/F’', fc = 1,2,..., whose heights are proportional 
to NE"^/kT. When added together, the peaks may give 
a highly irregular shape, as shown in Fig. 5. 


( 15 ) 


V. THE INTERFERENCE MECHANISM 


P'^(r) = Nw + E FjA, E^^j), (13) The task of evaluating P'^{t) is particularly straight- 

j forward, since we need not follow the evolution of the 
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p'^(x-ct-ct^) 


FIG. 6. (Colo r on line) Real part (dashed) and modulus (solid) 
of ipn in Eq.(16l transmitted across a barrier supporting a 
single metastable state, with p'^T/c = 0.05 and A(p^)T/c = 
0.11. Also shown by the dot-dashed line is the modulus of 'ipn 
as give by Eq.([T^ 


wave function in the barrier region, an only require the 
final overlap matrix Tmn- By filtering momenta of a ■i/'n, 
a narrow resonance at EJ produces a nearly monochro¬ 
matic transmitted state ibZn with E E^. The state is 
broad in the coordinate space, and the overlaps between 
and whose relative delay is (m —n)T, contains a 
factor exp[iE''j{m — n)T\. For several well-separated reso¬ 
nances, and with j ^ f have different energies, 
and are practically orthogonal, so that transmissions via 
different metastable states are mutually exclusive events. 
In this way, summation over all summation over all trans¬ 
mitted modes in Eq.([^ produces interference structures 
shown in Figs. 3-5. 

We illustrate this with a simple example, by considering 
non-spreading wave packet states with a linear dispersion 
law, E{p) = cp. Assuming the Breit-Wigner form for the 
transmission amplitude, T{p) = (F/[(£’ — E'') + IF], and 
putting A{p) « we have 

(x, t) » — ct — ct„) (16) 

c 

where 


^(y) = Si-y) exp{ip^y + Fy/c), 


(17) 


and d(y) = 1 for y > 0 and 0 otherwise. Thus, ■i/'n 
in Fig.6 has a sharp front followed by a long expo¬ 
nential tail, which allows to overlap, (V’mlV'n) 
exp[—iE^{tm — tn)], even if ipn didn’t, Imn = 0 for m n. 
With two, or more, resonances involved would con¬ 
tain several contributions of the form (161, one for each 
metastable state. 

Our analysis of a single particle, prepared in an exotic ini¬ 
tial state, helps us to gain an insight into the interference 
effects accompanying scattering of several identical parti¬ 
cles. Recently we considered [5] the case of two fermions 
or bosons, emitted in the same wave packet state, with 
a time delay r between the emissions. For a barrier with 
two resonance levels, the 2-particle transmission proba¬ 
bility P^(2,2), considered as a function of r, exhibited 


oscillations with a frequency Aw = E 2 — E^. The oscil¬ 
lations disappear if the particles can be distinguished. 
We note that in both problems we only require the initial 
and final overlap matrices, Imn, and Tmn [S]- They are, 
however, used differently. Whereas in the single particle 
case, considered here, the correction SP^{t) = 2ReTi2 
contains all barrier frequencies, EJ, the correction to 
P^(2,2) depends on |Ti 2 p, and oscillates only with Aw 
What makes the two cases similar, is that the particle, 
or particles, are distributed between wave packet modes 
tpi and ip 2 - For a single particle, this is readily seen 
from Eq.([^. A symmetrised (antisymmetrised) state of 
two uncorrelated particles, I 12 = I 21 = 0 is given by 
'i>{xi,X2,t) = [tpi{xi,t)i;2{x2,t)±'ijji{x2,t)tj;2{xi,t)]/V2, 
which also implies that particle 1 is simultaneously 
present in both ipi and ^’ 2 , albeit in a different manner. 
In both cases, the physical origin of an oscillatory pattern 
in the transmission probability is the overlap between dif¬ 
ferent '(/'n ) acquired in the barrier, and the phases carried 
by different ipn, launched at different times. 


VI. TRANSMISSION OF A MIXED CAT STATE 


Before concluding, we briefly discuss transmission of 
a mixed cat state with two components, tjji{x,t) and 
'Ip 2 {x,t), (V'mIV’m) = ^mn, rn,n = 1,2. The system is 
prepared in the following way: with a probability p/2 it 
is in one of the states ipi{x,t) and ip 2 {x,t), and with a 
probability (1 — p) it is in their coherent superposition, 
[ilji{x,t) + ip2{x,t)]/'/2. The incident density matrix is, 
therefore, given by 

P{X,X') = [-lpi{x,t)i:l{x',t) + 'llj2{x,t)'ip2ix' ,t)]/2 + 

(1 - p)[tpi{x',t)ilj 2 {x,t) +'il; 2 {x,t)'iljlix,t)]/ 2 . (18) 


For the transmission probability we have 

= [u;i + W2]/2 + (1 - p)i?e[Ti2(T)] (19) 


Thus, for a pure state, p = 0 we recover Eq.([^. As p 
increases, the last interference term in Eq.(19) becomes 
smaller, and finally vanishes for the incoherent combina¬ 
tion of the two states, p = 1, where P^ = (wi -I- W2)/2. 
This simple result is easily extended to the case where 
the initial state has three, or more, components, N > 2. 


VII. CONCLUSIONS AND DISCUSSION 

In summary, various components of the same one- 
particle state, well separated initially, may coincide inside 
a scatterer, provided the particle is detained there for an 
appreciable period of time. Then the interference, result¬ 
ing from this ’’pile up” effect, may significantly change 
the tunnelling probability The effect requires that 
transmitted modes are significantly broadened in the co¬ 
ordinate space, and is absent in tunnelling across a rect¬ 
angular barrier. It is present in resonance tunnelling. 
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where considered as a function of temporal delays 
between the components, oscillates with internal frequen¬ 
cies of the barrier. We have shown that the interference 
patterns predicted for the resonance transmission of sev¬ 
eral identical particles i. 0, have a similar origin, both 
resulting from the particle being distributed, in one way 
or another, between different wave packet components. 


The interference patterns are washed out if the initial 
state is mixed, rather than pure. The proposed type of 
’’interferometry in the time domain” is within capability 
of modern experimental techniques mm- 
Support of the Basque Government (Grant No. IT-472- 
10), of the Ministry of Science and Innovation of Spain 
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